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ABSTRACT 

We present the mass functions of actively accreting supermassive black holes over the redshift range 
0.3 < z < 5 for a well-defined, homogeneous sample of 15,180 quasars from the Sloan Digital Sky 
Survey Data Release 3 (SDSS DR3) within an effective area of 1644 deg^. This sample is the most 
uniform statistically significant subset available for the DR3 quasar sample. It was used for the DR3 
quasar luminosity function, presented by Richards et al., and is the only sample suitable for the 
determination of the SDSS quasar black hole mass function. The sample extends from j = 15 to i 
= 19.1 at z < 3 and to i = 20.2 for z > 3. The mass functions display a rise and fall in the space 
density distribution of active black holes at all epochs. Within the uncertainties the high-mass decline 
is consistent with a constant slope of /3 « —3.3 at all epochs. This slope is similar to the bright end 
slope of the luminosity function for epochs below z = 4. Our tests suggest that the down-turn toward 
lower mass values is due to incompleteness of the quasar sample with respect to black hole mass. 
Further details and analysis of these mass functions will be presented in forthcoming papers. 
Subject headings: cosmology: observations - galaxies: active - galaxies: mass function - quasars: 
emission lines - quasars: general - surveys 



1. INTRODUCTION 

We live in exciting times. Since their discovery, the 
power-house of quasars was believed to be actively ac- 
creting supermassive black holes (e.g., Salpeter 1964; 
Zeldovich & Novikov 1964). Owing to recent advances, 
not only have we been able to confirm the existence of 
such dense, powerful objects (e.g., Ghez et al. 2005; Gen- 
zel et al. 2003; Harms et al. 1994) the determination of 
their mass has become a common focus of many stud- 
ies, in spite of the typically non-trivial nature of this 
process. The reasoning is quite simple: (a) for local 
galaxies, the early studies revealed a close connection 
between the black hole mass and the mass, luminosity, 
and velocity dispersion of the spheroidal component of 
the host galaxy (e.g., Magorrian et al. 1998; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000a, 2000b) suggesting 
an intimate connection between their evolutions and (b) 
recent semi-analytical and hydrodynamical simulations 
(e.g., Granato et al. 2004; Springel et al. 2005; Catta- 
neo et al. 2006; Kang et al. 2006) show that black hole 
evolution and activity may have an important influence 
on the properties of massive elliptical galaxies. Mapping 
the demography and mass of the black holes in the cen- 
ters of galaxies across the history of the Universe is a 
first step to understanding this connection. Since a frac- 
tion of the matter accreted onto black holes as they grow 
is converted into radiation, the cosmic evolution of the 
quasar luminosity function helps constrain the accretion 
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and growth history of black holes and can therefore also 
help us understand the connection with galaxy evolu- 
tion. Yet certain degeneracies limit what we can learn 
from the luminosity functions (e.g., Wyithe & Padman- 
abhan 2006) owing to the unknown degree of the radia- 
tive efficiency and the mass accretion rate and how these 
parameters evolve and depend on black hole properties. 

In principle, these degeneracies can be broken by com- 
bining the luminosity and mass functions. Recently, 
Richards et al. (2006, hereafter R06) presented the lu- 
minosity function and its cosmic evolution from redshift 
5 to redshift 0.3 for a well-defined, homogeneous sub- 
set of 15,180 quasars from the Sloan Digital Sky Sur- 
vey (SDSS) Data Release 3 (DR3) quasar catalog of over 
46,400 quasars for which the quasar selection function 
is well-determined. Therefore, the determination of the 
mass function for this specific quasar sample is of high 
interest. In this Letter we present this mass function of 
actively accreting black holes and its cosmic evolution 
from redshift 5 to 0.3. 

To determine the mass of the black holes we adopt 
the scaling relations based on the broad emission line 
widths and nuclear continuum luminosities (e.g., Vester- 
gaard 2002; Warner et al. 2003; Vestergaard & Peterson 
2006) because (a) they are anchored in robust black hole 
mass determinations of active black holes in the nearby 
Universe (Peterson et al. 2004; Onken et al. 2004), (b) 
several lines of evidence exist in favor of our application 
of this method to the most distant active black holes 
(see Vestergaard 2007 for details), and (c) the associated 
uncertainties on the mass estimates are among the low- 
est for the available mass estimation methods for distant 
sources. A cosmology of Hq = 70 km s^^Mpc^^, VIa = 
0.7, and = 0.3 is used throughout. 

2. SPECTRAL MEASUREMENTS 

As noted above, we use the quasar sample presented by 
R06 because it is well-defined and has a well-understood 
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selection function; this is not the case for the fuU DR3 or 
the DR5 sample. Therefore it is currently also the only 
sample suitable for the determination of the black hole 
mass function. The reader is referred to R06 for details 
on the sample and its selection. 

In order to estimate the mass of the central black hole 
in each quasar we need to measure the widths of each 
of the H/3, Mgii, and the Civ emission lines and the 
monochromatic nuclear continuum luminosity near these 
emission lines. Since quasar spectra contain contribu- 
tions from other emission components which often 'con- 
taminate' the line and continuum components of interest, 
we model each of these emission contributions so to ob- 
tain reliable measurements following a decomposition of 
the modeled components. In the following we provide 
only a very brief summary of our data handling. A more 
detailed account of our processing of the SDSS spectra 
will be presented by M. Vestergaard et al., (in prepara- 
tion). 

Each spectrum is corrected for Galactic reddening 
based on the E(_B — V) value relevant for the quasar 
(Schneider et al. 2005) and the Galactic extinction maps 
by Schlegel et al. (1998). The continuum components 
were modeled using: (a) a nuclear power-law continuum, 
(b) an optical-UV iron line spectrum (Veron et al. 2004; 
Vestergaard & Wilkes 2001), (c) a Balmer continuum, 
and (d) a host galaxy spectrum (for objects at z < 0.5 
for which the galaxy contribution is strongest and is best 
characterized; Bruzual & Chariot 2003). 

The continuum components were subtracted and the 
emission lines were then modeled with multiple Gaus- 
sian functions so to obtain smooth representations of the 
data. This allows us to eliminate most narrow absorption 
lines, noise spikes, and the contribution from the Narrow 
Emission Line Region, which is by far the strongest in 
the optical region. A single Gaussian component, not to 
exceed a FWHM of 600 kms~^, was used to model and 
subtract the narrow line components. The line widths of 
up, Mgll, and Civ were measured for all emission lines 
with model line peaks above three times the median flux 
density error across the emission line (i.e., > 3ct peaks). 

All Mg II and C iv profiles with strong absorption were 
discarded from further analysis. They were identified by 
visual inspection of the spectra of the quasars listed in 
the catalog of Trump et al. (2006) to have absorption 
and of the quasars with redshifts between 1.4, when C iv 
enters the observing window, and 1.7. These quasars are 
prone to have strong C IV absorption but this redshift 
range is not covered by Trump et al. (2006). 

Of the 15,180 quasars on which the DR3 quasar lumi- 
nosity function (R06) is based, black hole mass estimates 
were possible for 14,434 quasars (95%). 

3. BLACK HOLE MASS ESTIMATES 

The specific mass scaling relations we use are equa- 
tions 5 and 7 of Vestergaard & Peterson (2006) based on 
the FWHM of H/3 and C iv and a new relationship for 
the Mgii emission line based on high-quality data from 
the SDSS DR3 quasar sample (M. Vestergaard et al., in 
preparation) . For each quasar a black hole mass estimate 
is computed for each of the H/3, Mgii, and C iv emission 
lines rendered suitable for mass estimates (§[2|). The final 
black hole mass was computed as the variance weighted 
average of the individual mass estimates. The mass esti- 
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Fig. 1. — Redshift distribution of the black hole masses of our 
quasar sample selected from SDSS DR3. Of the 15,343 quasars 
detected in the selected sky area at 2 > 0.2, mass estimates were 
possible for 14,584 quasars. Of these, 14,434 quasars have redshifts 
between 0.3 and 5.0. The median mass in each redshift bin shown 
in Figure[2]is marked with the (red) open box. The median propa- 
gated black hole mass measurement error is shown in the lower part 
of the diagram. For reference, the (green) dashed curve indicates 
the faint SDSS flux limits folded with a line width of 2000 kms'^ 

mates based on the individual emission lines match well 
within their errors. 

Figure [1] shows the redshift distribution of the black 
hole masses. Its detailed shape is mainly a consequence 
of the survey properties and limitations: (a) the scarcity 
of quasars at 2.2 < z < 3.0 is due to inefficient quasar se- 
lection, discussed later and (b) the lower boundary is de- 
termined mainly by the faint survey flux limits {i < 19.1 
mag at z < 3.0; i < 20.2 mag at z > 3). While the bright 
flux limit (i > 15 mag) of SDSS is expected to eliminate 
the most massive black holes at z < 1, this does not have 
a strong impact on the current study. Judging from the 
comparison made (Jester et al. 2005) between the Bright 
Quasar Survey (hereafter BQS, Schmidt and Green 1983) 
and SDSS the surface density of quasars brighter than i 
~ 15.0 is such that at least 6 such bright quasars are ex- 
pected in the 1622 square degrees of the current sample; 
this is a lower limit since the BQS itself is incomplete. If 
these bright quasars follow the same redshift distribution 
as the BQS quasars detected by SDSS, they will mostly 
lie at z < 0.25, below the minimum redshift bound (0.3) 
of the mass functions. Also, no upper limit in line width 
is imposed that would place an upper bound on the black 
hole masses. While we do deselect low quality emission 
line profiles (see § ^ this process applies to all profile 
widths (and thus all black hole masses). Given our large 
sample size (nearly 15,000 quasars) and that we know 
of no selection effects that would systematically deselect 
the most massive black holes, it is fair to conclude that 
the upper bound in black hole mass is real. The mass 
distribution in Figure [T] therefore confirms with larger 
numbers the results from previous studies (e.g., Corbett 
et al. 2003; Warner et al. 2003; Netzer 2003; McLure 
& Dunlop 2004; Vestergaard 2004; KoUmeier et al. 2006; 
Netzer & Trakhtenbroot 2007; see also Vestergaard 2006, 
2007; Shen et al. 2007) that find the detected quasars at 
high redshift have black hole masses between a few times 
lO^M© to 10^° Af© with typical mass of order 10^ Mq. 

4. BLACK HOLE MASS FUNCTIONS 

The differential quasar black hole mass function, 
5'(A/bh,z) is the space density of quasar black holes (i.e., 
the number of black holes per unit comoving volume) per 
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Fig. 2. — Mass functions of active supermassive black holes in the 
quasar sample at a range of redshifts (open diamonds; error bars 
denote the Poisson statistical uncertainty). The mean redshift of 
each redshift bin is labeled. The boundaries of the redshift bins 
are: 0.3, 0.68, 1.06, 1.44, 1.82, 2.2, 2.6, 3.0, 3.5, 4.0, 4.5, 5.0. In 
each panel the z = 0.49 and 2.01 mass functions are also shown as 
(black) dotted and (red) dashed curves for reference. The vertical 
(green) dashed and (blue) dotted lines mark the masses of two 
fiducial black holes at the mean redshift of the z-bin, emitting a 
luminosity equal to the survey flux limit, and around which the 
broad emission line gas is moving at 2000 kms""^ and 3500 kms~^ 
(FWHM), respectively. Incompleteness is expected for black hole 
masses below the 2000 kms~^ equivalent fiducial mass value, but 
may also affect bins of higher mass values. See text for details and 
discussion. 

unit black hole mass as a function of black hole mass and 
redshift. We calculate this space density in different mass 
and redshift bins using the 1/Vq method following War- 
ren, Hewett, and Osmer (1994), where Va is the accessible 
volume, defined by Avni and Bahcall (1980). We adopt 
the method of Fan et al. (2001) of including the selec- 
tion function of the quasars: we interpolate the selection 
function between the grid points defined by the quasar 
luminosity and redshift and integrate over the selection 
function in determining Va\ see equation (6) of Fan et al. 
The SDSS DR3 quasar selection function for our quasar 
sample is published by R06. The mass function (and its 
statistical uncertainty) is given by 



v['(<M>,<2>) = ^— -i 



E 



1 



-,1/2 



(1) 



(2) 



where < M > and <z> are the average mass and red- 
shift of each mass bin of size AMbh and redshift bin of 
size Az and index i refers to each individual quasar in 
this mass and redshift bin. In Figure[2]we show the mass 
functions at a range of redshifts. For better comparison, 
the redshift bins are chosen to be those used for the lumi- 
nosity function for this sample (R06). The data points 
with large errors are typically due to very few quasars 
residing in those mass bins. 

5. DISCUSSION AND SUMMARY 

For each redshift bin in Figure [5] the mass function 
displays a rise and fall in the mass distribution. The 
decline toward higher masses is expected to be intrinsic 



to the quasar population (§[3|); while the uncertainty in 
the mass estimates will affect the accuracy of the slope 
value, it is too small to be responsible for the decline. We 
determine^ a typical high-end slope (3 of about —3.3 (^' oc 
M^) with uncertainties between 0.4 to 1.3 for the mass 
functions below a redshift of 4. At the mean redshifts of 
4.25 and 4.75 the high-end slope 'flattens' to /? « -2.9 ± 
1.1 and f3 w —1.9 ± 1.7, respectively. The mass function 
at redshift 4.75 is based on very few black holes at each 
mass bin, especially for Mbh > 10^'^ Mq, so the shape 
is somewhat uncertain. The mass functions are clearly 
consistent with a constant slope across all epochs. This 
is also easily verified in Figure [2] by the reference curves 
for the mass functions at 2=0.49 (black dotted) and at 
z= 2.01 (red dashed) which also show the space density 
of high mass black holes slightly increases at 2 > 1 . 

The slope of the mass functions is consistent with 
the bright end slope of the quasar luminosity function 
(/3 w —3.3; Croom et al. 2004; R06) to within the uncer- 
tainties with the exception that the latter flattens above 
a redshift of 4 (to /3 = -2.5; Schmidt et al. 1995; Fan 
et al. 2001). Another difference is the relative normal- 
izations. The mass functions exhibit a change in space 
density for a given mass of only a factor of a few between 
2 = 0.5 and 2 = 2 (cf. the dotted and dashed curves in 
Fig [2]), while the space density of quasars brighter than 
~ —27 mag drops 2 orders of magnitude or more between 
redshifts 2 and 0.5 (Fig. 18, R06). While there is no 
one-to-one correspondence between a given mass value 
and a given luminosity value since black holes radiate 
at a range of (Eddington) rates, the different space den- 
sity changes show that the cumulative volume density of 
L/M (oc L/LEdd) decreases toward low redshift. This 
confirms what we already know based on the declining 
space density of active nuclei at later epochs (e.g., R06): 
low redshift black holes are typically less active. 

In Figure [21 a characteristic peak shift to higher mass 
values with redshift is seen. It is important to establish 
the reality of these peaks and their shifts, since this would 
have intriguing cosmological implications. For one, it 
would be a clear manifestation of cosmic down-sizing: 
the most massive black holes were more actively accret- 
ing early in the universe and with time progressively less 
massive black holes dominate the population of actively 
accreting black holes. Alternatively, the down-turn to- 
ward lower mass values may simply be due to incom- 
pleteness of the sample, as the source selection is based 
on the source luminosity and broad-band colors, not on 
black hole mass. As a result, we will not obtain a sharp 
lower boundary in mass (as we do in the luminosity distri- 
bution) but a decrease in space density below the survey 
flux limits, owing to the distribution of black hole masses 
at a given luminosity. Unfortunately, the SDSS does not 
go deep enough for us to make meaningful cuts in black 
hole mass for the current sample, as will be clear later. 

The best way to assess which part of the down-turn 

^ The slope of the mass functions at 2; < 2 is determined for 
masses above 9.2 dex and at z > 2 for Mbh > 9.6 dex, since these 
mass bins are deemed (essentially) completely populated. To ac- 
count somewhat for the uncertainties in the Mbh values and that 
the mass bin population varies, the uncertainty in a given mass bin 
was estimated as the propagated uncertainty in the absolute mass 
estimates (e.g.. Table 5 of Vestergaard & Peterson 2006) given the 
Poisson population statistics. 
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is intrinsic requires detailed simulations that include the 
uncertainties in the mass estimates, the selection func- 
tion, and the flux limits. Since this is non-trivial we 
will address this in a forthcoming paper (B. Kelly et al., 
in preparation). Here, we instead perform crude checks 
as follows. Firstly, in Figure [2] we mark the masses of 
a fiducial source at the mean bin redshift with a con- 
tinuum luminosity equivalent to the flux limit at that 
redshift and a broad line width of 2000 kms^^ (green 
dashed vertical line) and 3500 kms^^ (blue dotted verti- 
cal line), respectively; 2000 kms^^ is the canonical lower 
line width for which a source is a bona-fide quasar and 
3500 kms"'^ is the median linewidth for nearby quasars. 
There are objects below the (green) dashed line since the 
observed i-band magnitude contains emission contribu- 
tion in addition to the power-law continuum luminosity 
used for the mass estimates (see §[5]). Nonetheless, these 
fiducial masses do yield a guideline location of the SDSS 
flux limits in the mass parameter space, and their loca- 
tion close to the peak of the mass distribution suggests 
that the peak shift is simply due to the survey flux limits. 

Given the non-trivial relationship between the z-band 
magnitude and the monochromatic luminosity, we also 
performed a simple simulation with the aim of mini- 
mizing selection effects. We generated a mock cata- 
log of quasars covering redshifts between 0.2 to 5.0 and 
bolometric luminosities between 10^^ ergs~^ and fO*** 
ergs~^. For each (L, z) grid point we computed the dis- 
tribution of black hole masses by applying the typical 
observed line width distribution for these L and z val- 
ues. The mass functions based on this mock catalog (not 
shown) do not display the narrowly peaked distributions 
seen in Figure[2l but display a slower continued rise in the 
distribution below the high mass fast decline and a sharp 
downturn at the lowest masses which coincides with the 
lower luminosity bound. No downturn is seen at the low 
mass end if a sharp cut in black hole mass is imposed. 
Hence, this mock catalog confirms the earlier indication 
that the low-mass downturn in the observed mass func- 
tions is a consequence of the relatively narrow luminosity 
range probed by SDSS broadened by the line width dis- 
tribution, especially at higher redshift. We thus conclude 
that the peaks and their shifts with redshift are mainly 
due to incompleteness of active black holes near the sur- 
vey flux limits. These crude checks are, however, unable 



to verify which part of the turnover is trustworthy. The 
mass functions in the lower redshift bins display a slower 
turnover, part of which may well be real. B. Kelly et al., 
(in preparation) will present more advanced simulations 
that address these and related issues. 

The space density decline at z « 2.8 relative to z = 2.4 
and z = 3.25 in Figure [2] is caused by a much lower selec- 
tion efficiency (by a factor of 8; Fig. 6 of R06) as quasars 
at these epochs coincide with the stellar loci in the se- 
lected color-spaces. This suggests that the completeness 
level is slightly overestimated at these epochs. This sit- 
uation is unfortunate since the space density of quasars 
is known to peak between redshifts of 2 and 3 (e.g., Os- 
mer 1982; Warren et al. 1994; Schmidt, Schneider, & 
Gunn 1995; Fan et al. 2001). It is thus not possible with 
the current sample to get a complete picture of how the 
space density distribution of black holes behave during 
this important epoch. The Large Bright Quasar Sur- 
vey (hereafter LBQS; Hewett et al. 1995) of about 1000 
quasars, extending to a redshift of 3, offers an oppor- 
tunity to address these issues further, since this sample 
is not selected based on broad-band color. The LBQS 
quasar mass function will be presented by Vestergaard & 
Osmer (in preparation) and discussed in concert with the 
DR3 mass function by M. Vestergaard et al. (in prepa- 
ration), where the mass functions will also be compared 
to previous work. 
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